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We investigate experimentally a deterministic underdamped Josephson vortex ratchet - a fluxon- 
particle moving along a Josephson junction in an asymmetric periodic potential. By applying a 
sinusoidal driving current one can compel the vortex to move in a certain direction, producing 
average dc voltage across the junction. Being in such a rectification regime we also load the ratchet, 
i.e., apply an additional dc bias current /dc (counterforce) which tilts the potential so that the 
fluxon climbs uphill due to the ratchet effect. The value of the bias current at which the fluxon 
stops climbing up defines the strength of the ratchet effect and is determined experimentally. This 
allows us to estimate the loading capability of the ratchet, the output power and efficiency. For 
the quasi-static regime we present a simple model which delivers simple analytic expressions for the 
above mentioned figures of merit. 

PACS numbers: 05.40.-a, 05.45.-a, 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

The discovery of Brownian motion gave birth to the 
idea of extracting useful work out of random motion. 
As Richard Feynman et al. demonstrated draw- 
ing energy from equilibrium thermal fluctuations (white 
noise) is forbidden by the second law of thermodynamics. 
The extraction of work out of non-equilibrium or time- 
correlated noise (colored noise) is possible in the so-called 
ratchet systems. Such systems incorporate an asymmet- 
ric periodic potential and have been in the focus of in- 
vestigation during the last two decades in various imple- 
mentations. Ratchets based on Josephson junctions, in 
particular exploiting the motion of the Josephson phase 
in SQUIDs or vortices in long Josephson junctions (LJJs) 
have been suggested and tested experimentally . 

These systems have some advantages over other sys- 
tems: (I) directed motion results in an average dc volt- 
age which makes ratchet operation easily accessible in ex- 
periment; (II) Josephson junctions arc very fast devices, 
i.e., they can be operated in a broad frequency range 
from dc up to 100 GHz which allows them to capture a 
lot of spectral energy; (III) both underdamped and over- 
damped systems can be investigated by proper junction 
design and the variation of the bath temperature. 

The deterministic underdamped Josephson vortex 
ratchet (JVR) in which a Josephson vortex (fluxon) 
moves along a LJJ was implemented by our group ear- 
lier [l5j . A strongly asymmetric tunable potential was 
created using a current injector [ll|, IH| • The periodic- 
ity of the potential is provided by the annular geome- 
try of the LJj [l7l - [20j . The fluxon was injected into the 
annular LJJ also in a controllable way using a pair of 
tiny current injectors pll . [22j . The directional motion 
of a fluxon was detected by measuring the (averaged) dc 
voltage across the junction which is, due to the Joseph- 
son relation, proportional to the average velocity of a 
fluxon and was reaching values as high as 0.9cq, where 



Co is the Swihart velocity (maximum possible velocity of 
a fluxon). However, in these experiments the JVR was 
idle, i.e. was not delivering any rectified power to a load. 

In this paper we investigate several figures of merit of 
such a JVR relevant for applications: rectification win- 
dow, maximum dc counterforce against which the ratchet 
can still work, output power and efficiency in both quasi- 
static and non-adiabatic regimes. 

The paper is organized as follows: In Sec. |TT] we de- 
scribe the equations for the dynamics of the Josephson 
phase in a LJJ with a ratchet potential created by an 
injector current. The sample design and characterization 
are presented in Sec. IIIII Experimental and analytical 
results for quasi-static driving frequencies are presented 
in Sec. II VI The Sec. IV1 covers experimental and numeri- 
cal results for the non-adiabatic drive. Sec. I VII concludes 
this work. 



II. THEORY 

Our system consists of an annular long Josephson junc- 
tion (ALJJ) equipped with injectors to create an asym- 
metric potential and to insert a fluxon, see Fig. [TJ The 
dynamics of the Josephson phase in the system can be de- 
scribed by the following perturbed sine-Gordon equation 

<Pxx - 4>tt - sin<j) = a^t-J- Tmj(x) - £(*), (1) 

where (j)(x, t) is the Josephson phase and subscripts x and 
t denote derivatives with respect of space and time, re- 
spectively. The curvilinear coordinate x along the LJJ is 
normalized to the Josephson penetration depth Aj and 
the time t is normalized to the inverse of the plasma 
frequency UJZ . The quantity a is the dimensionless 
damping parameter. 7 = j/j c , y m j(x) = j- m j(x)/j c and 
£(t) = j(t)/j c are the dc bias current density, injector 
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FIG. 1: (Color online) (a) Optical image of the ALJJ used in 
our experiments. Only one injector of the right pair is used, 
(b) Sketch of the ALJJ with a single injector 1 and a double 
injector 2. Bias leads for top and bottom electrode are not 
shown for clarity. 



current density and ac driving current density, respec- 
tively, all normalized to the critical current density j c of 
the LJJ. 7i n j(x) has zero spatial average and is used to 
create an asymmetric potential, 7 is used to apply an 
additional dc bias current while the device is operated in 
the rectification regime (see below). £(i) is a spatially ho- 
mogeneous deterministic (or stochastic) drive with zero 
time average. The ultimate aim of ratchet operation is 
to rectify to produce a nonzero voltage V oc 4> t 7^ 0, 
which is independent on x. 

In the absence of the right-hand side, the solitonic so- 
lution of Eq. fjj) is a Joscphson vortex (sine-Gordon kink) 
(f>(x) = 4arctan [exp {[x — xo(t)} j\J\ — u 2 }] situated at 
x and moving with velocity u = dxQ(t)/dt [23[. The 
right-hand side of Eq. fJJ) is usually considered as a per- 
turbation |24[. It does not change drastically the vor- 
tex shape, but defines its dynamics, e.g., its equilibrium 
velocity [24| . Such an approximation essentially treats 



the vortex as a rigid object, and its dynamics can be 
reduced to the dynamics of a relativistic undcrdarnped 
point-like particle [l2j (cf. the nonrelativistic case f25j). 
In these terms, the ratchet should rectify £(i) to produce 
a nonzero average velocity «^0. 

We implement the asymmetric potential using a sin- 
gle injector (further called injl), see Fig. [TJ A suit- 
able current profile 7i n j (x) is equivalent to a non-uniform 
magnetic field h{x) along the LJJ, such that 7^(2;) = 
—h x (x). The normalized potential U(xq) felt by a fluxon- 
particle is given by U(xq) ~ — 2-7T (wj/Aj) h(xo), where 
Wj is the junction width [111 . 
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TABLE I: Parameters of the used junctions, r is the junc- 
tion radius, j c is the critical current density at 4.2 K, Aj is 
the Josephson penetration depth, I is the normalized junc- 
tion length. Ax and Aw describe the injector seperation and 
width (see text). 



Such a profile can be realized by applying a current 7; n j 
through a single injector of width Aw situated at x = 
Xi n ji and extracting the current over the same electrode 
along the rest of the LJJ, i.e., 



7inj(x) 



71 for \x - Xi n ji| < Aw/(2Xj) 

72 otherwise 



(2) 



such that 71 Aw + 72 (/ — Aw) = with the ALJJ cir- 
cumference L = 2irR, where R is the mean radius of 
the ALJJ. In this case, U(xq) looks like an asymmetric 
saw-tooth potential with the steep slope oc 71 and the 
gentle slope cx 72, i.e., the asymmetry depends on the 
width Aw compared to the junction length L, see Fig. [1] 
The amplitude of the potential can be varied by chang- 
ing Jjnji = ^YiAwwj, which, in principle, also allows the 
operation as a flashing ratchet. Here, we focus on the 
rocking ratchet only, i.e., the potential oc 7i n ji is (al- 
most) constant and £(t) 7^ 0. To apply the current Ji n ji 
we use one of the two injectors visible in Fig. [TJ The in- 
jector is attached to the bottom electrode. The current 
Jinji is injected via injl into the bottom electrode of the 
ALJJ and is used to create the asymmetric potential as 
described above. 

The periodicity of the potential is provided by using an 
annular LJJ (ALJJ) and (7i n j(x)), = 0. Note, that the 
latter condition is automatically satisfied because Ji n ji 
is applied to the same electrode so that jij c Awwj = 

l23c{L - Aw)wj = /injl- 

In addition, our JVR has a pair of current injectors 
(further called inj2) separated by a distance Ax and at- 
tached to the top superconducting electrode, see Fig. [TJ 
They are used to insert a fluxon (Josephson vortex) in 
the ALJJ fHHl. 



III. SAMPLES 

We investigated several Nb/Al-AlO^/Nb junctions 
with different parameters. Here, we report the results 
obtained for two samples with parameters summarized 
in Table [TJ For all samples, Aj 3> wj = 5 /im and 
Aj ^> Aw = Ax (see Table [TJ), i.e., we can treat our 
ALJJs as one-dimensional and inj2 as an ideal discon- 
tinuity [22| . The maximum revolution frequency for a 
fluxon inside an ALJJ equals = c /L with the Swi- 
hart velocity cq. The corresponding voltage (voltage of 
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FIG. 2: /c(^inji) dependence for sample C3. The range of 
/inji relevant for ratchet operation, i.e., the region where the 
picture of a fluxon-particle moving in a potential is valid, is 
shaded. 



the first fluxon step) is given as 



(3) 



with the normalized length I = L/Xj. 

All measurements have been performed at the temper- 
ature T = 4.2 K. All junctions showed good I-V char- 
acteristics (IVC) and symmetric 7 C (77) dependences (not 
shown). Their / c (iinj2) dependences look Fraunhofer-like 
in agreement with the theory [22I ] . The first minimum of 
this dependence corresponds to the phase being twisted 
by ±27r in a tiny region between the injector pair inj2 and 
to a free (anti)fluxon being inserted into the ALJJ region 
outside inj2. Thus, we insert a fluxon into the junction 
by choosing the corresponding value of I m j2- 

To calibrate injl, we measure 7 c (7; n ji) (see Fig. [2]). By 
measuring 7 c (7i n ji) without a fluxon inside the junction 
(-finj2 = 0) and comparing it to 7 c (/i n ji) measured with 
a fluxon inside the junction (ii n j2 oc A(f> = ±27r), the 
asymmetry of the ratchet is determined and the range of 
magnetic field (range of the amplitude of the ratchet po- 
tential) relevant to the ratchet operation is found (TlT[l5j . 
see Fig. [2] Inside this working range, the depinning cur- 
rent 7 C scales almost linearly with imji, but it is asymmet- 
ric for positive and negative direction of the bias current 
7 (driving force). Note, that like in Ref. [HI a residual 
pinning of the fluxon due to the finite inj2 sizes Ax and 
Aw is still visible in Fig. [5] [22j ■ 

Measuring an IVC with a fluxon inside the junction 
and 7j n ji / 0, i.e., an applied ratchet potential, a fluxon 
step appears on the IVC (not shown) corresponding to 
the rotation of a fluxon around the ALJJ with u ~ co. 
The depinning currents 7 C and return currents I r of the 
fluxon step depend on the polarity of the applied bias as 
well as on L 



injl- 



IV. QUASI-STATIC DRIVE 

For our experiments in the quasi-static 
regime, we apply a periodic bias current 
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FIG. 3: (a) Typical rectification curves V(/ a c) for Jmji = 
±200 /iA, ±400 fj,A and ±600 fiA (open, black and grey sym- 
bols, respectively), (b) Figures of merit for different potential 
amplitudes normalized to Iq and Vn, respectively. A/ rcc t is 
shown for both positive and negative potential amplitude. 



I(t) = t;(t)j c Lwj = I ac sin(27rt4) with the frequency 
v = 100 Hz <C and measure the rectification curve 
V^/ac) by averaging the voltage over 10 ms (1000 data 
points sampled at 100 kHz) — one period of the ac drive. 
For the junction E3 the dependence V(I ac ) = Vi ■ u(I ac ) 
is shown in Fig. [3] (a) for different values of /mjij 
for different amplitudes of the potential. Note, that 
the absolute maximum V , that can be obtained with 
a rectangular-shaped driving force, is given by 0.5 • V\, 
because the fluxon is driven only half a period in one 
direction. For a sine-like driving force, as in our case, 
the maximum voltage is somewhat smaller 11 [ . 



All V (lac) curves show similar features. For small I ac 
the driving force acting on a fluxon is not sufficient to 
push the fluxon out of the potential well in either direc- 
tion so that u oc V = 0. At . 7 ac > 7™" the bias is able 
to push the fluxon in one direction but not in the other, 
which results in u oc V / 0. At 7 ac > 7™^ t x , the junction 
switches into the resistive state generating a high positive 
or negative dc voltage. The latter regime is not discussed 
here as it has nothing to do with the JVR operation. 
Rarely, we also observed the typical ratchet behavior - 
a decrease of V at 7 ac > 7™^ t x when the driving force is 
able to overcome the potential barrier in both directions 
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and rectification, thus, drops significantly. We rarely ob- 
serve this regime since the asymmetry is so large that a 
negative fluxon step does not appear in most cases. 

For not very large potential heights oc -Zinji, i-e., when 
the perturbation theory is applicable, the values /™™ 
and /™cf grow approximately linearly with imji, see 
Fig. [UJb). Therefore, the size of the rectification window 
Afreet = I^t ~ ^rect a ^ s0 g rows approximately linearly 
with /inji- Further, V max grows with Imji, but it reaches 
its maximum value already for medium values of iinjij 
see Fig. [3Jb). At high potential amplitudes the recti- 
fication window is becoming smaller, probably because 
the ratchet potential cannot be considered anymore as a 
perturbation. 

The dependence shown in Fig. (3fb) is qualitatively 
the same for all measured samples and for both signs 
of the potential, as long as the junction length I is large 
enough. For I < 4 the asymmetry of the potential is 
strongly reduced because the potential is a convolution 
of the magnetic field profile created by injl and the fluxon 
shape pd|. The main result of Fig. E^b) is that in order 
to obtain large Vmax and large A/ rec t, one should operate 
the ratchet at large amplitudes of the potential, i.e., at 
large values of imp- 

In the quasi-static regime, all information about the 
ratchet operation can be derived from its IVC. Therefore, 
we can derive all figures of merit that we are interested in 
by using a model step- like IVC shown in Fig.[4|a), which 
is very similar to a real one. In this IVC that we express 
as 



V mod (I) = V 1 



-1 
0, 

+1, 



Ic- < 

Ic+ < 



< Ic+ 



(4) 



the fluxon dcpinning currents I c+ > and / c _ < have 
different values reflecting the asymmetry of the potential. 
The fluxon step is roughly approximated by a vertical 
step with infinite height. An applied current I(t) = I c \ c + 
I ac sin (ut) consists of a dc current and an ac current with 
frequency uj. 

To obtain a rectification curve V(Iac) we integrate the 
instant voltage over one period T = 2it/uj of ac drive, 
i.e., 



V(hc) = i / 



(5) 



Within our simple model, this can be integrated ex- 
plicitcly resulting in 



V = 
where 
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FIG. 4: (Color online) (a) Model step-function-like IVC © 
with I c + = O.l/o and I c - = — 0.4Jo, Io = j c wjL is the intrin- 
sic critical current, (b) V (lac) curves calculated using Eq. ([6]). 
(c) Pin(Tac) dHJ. (d) Pout (iac) CC2> ■ (e) power efficiency 
j] = -Pout/Pin for different values of i dc = 0, 0.05, 0.10, 0.14. 
(f) istop(iac), see Eqs. and (|12|l with saturation value 

/stop = -0.15 (0. 
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describe the part of V rectified during the positive and 
the negative semi-period of ac drive, respectively. 

The V(I ac ) calculated in this way is shown in Fig.rjjb). 
First, for 7 dc = our simple model gives 7™t = P+' 
^rect x = —Ic- and A/ rect = —I c - — I c +- For the chosen 
values of |7 C _| > |/ c +| the rectified voltage V > 0. 

Second, we apply 7 dc to try to stop the ratchet op- 
eration at a given value of 7 ac . The sign of 7 d c should 
be opposite to the sign of V, i.e. idc < in our case. 
Rectification curves V(I ac ) calculated using Eq. © at 
different values of Tic < are also shown in Fig. 0|b). 
One can see that the rectification window shrinks, i.e., 
close to the edges of the original window the ratchet is 
not strong enough to work against Idc- With applied Idc 



rin in 

rect 


= Ic+- 


Idc 


max 
rect 


= Idc — 


Ic- 




= 2/dc - 


Ic- 



Ic+, 



(8) 
(9) 
(10) 



which can also be seen from Eq. ©. 

One can also take a different point of view and study 
how strong the ratchet is at each particular 7 ac , i.e., 
one varies Idc at fixed 7 ac . The value of Idc at which 
the ratchet stops moving or starts moving backwards 
(V = or changes sign) is called stopping force (or 
stopping current) 7 s top. The dependence 7 s top (7 ac ) is 
shown in Fig. @|f ). One can see, that naturally 7 s top = 



up to the driving amplitude 7 a 



because the 



ratchet didn't even start working yet, so no force is 
needed to stop it. Then Tstop grows linearly to a value 
of i™« = A/ rcct /2 = (-J c _ - J c+ )/2, which is reached 
at I ac = (— J c _ + 7 c +)/2, i.e., at the middle of the rec- 
tification window. In this regime, the ratchet rectifies 
because during some fraction of the positive semi-period 
I(t) exceeds I c + and never exceeds 7 C _, see Fig. |U[a). 
An additional Idc "shifts" the origin of ac oscillations 
down along the /-axis. Thus, for Idc < |7 ac | — I c +, the 
ac sweep does not exceed the depinning current I c + and 
rectification vanishes. Therefore 



I 



stop 



= I he I - I c + for |7 ac | > I, 



c+, 



(11) 



which corresponds to the linear part in Fig.fJJf). The be- 
havior described by Eq. (|11[) does not hold for arbitrarily 
large I ac . When |/ ac | > (7 C+ — 7 c _)/2, at some amplitude 
of Idc the ac sweep also exceeds the negative depinning 
current 7 C _, which results in a drastic decrease of rec- 
tification and also 7™^ t x given by Eq. (JSJ) , see Fig. HJb) . 
At 7dc = (Ic- + 7 c +)/2, the ac sweep sees symmetric 
depinning currents and rectification vanishes. Thus 



7 



stop 



(7 C _ + 7 c+ )/2 for |7 ac | > (7 C+ - 7 c _)/2, (12) 



corresponding to a saturation of /stop in Fig. fjjf) 
and a completely closed rectification window on the 
^(7 a c,7stop) plot. For 7 dc < (7 G _ + 7 c+ )/2, the rectified 
voltage V changes sign. 

Let us now discuss a power balance in our ratchet. 
Since 7d c = const., the average dc output power 



1 

T 



V(t)I dc dt = VI dc < 0, 



(13) 



where V(I ac ,Idc) is given by Eq. ([6]). P out is negative 
because the ratchet delivers the power to the dc source 
(load) instead of consuming it. Simultaneously, the input 
power is given by 



P, 



V(t)I ac sm(ut) dt > 0. 



(14) 



Within our simple model, this can be integrated ex- 
plicitcly resulting in 



P, 



where 
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(15) 

(16a) 
(16b) 



describe the part of P n consumed during the positive and 
the negative semi-period of ac drive, respectively. 

Both Pn(7 a c) and P ou t(7 a c) at different values of 7d c are 
shown in Fig.UJc) and (d). One can see that Pi n (7 a c) has 
two characteristic branches corresponding to the dissipa- 
tion during the positive semi-period and during both the 
positive and negative semi-periods. In contrast, P ou t(7 ac ) 
is non monotonous and has an extremum given by 



y[C-(7 dc )]7 dc 



— arcsm 

7T 



f7 dc 
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\Idc 





(17) 
non- 



at 7 ac = 7™ c a t x (7 dc ) ©, which depends on 7 du 
monotonously. One can see that P ou t = not only at 
7 dc = 0, but also at 7 dc = 7 st o P where V = 0, c.f. Eq. (p~3|) . 
Therefore, the maximum power is reached for some in- 
termediate values of 7 st0 p < 7d c < 0. The exact value 
can be derived by looking for the extremum of expres- 
sion ([T7|) with respect to 7d c . It is reached for 7?^ which 
is a solution of the following transcendental equation 



7 C + — 7 C 
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i • MdT 
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= 0. 



(18) 
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For our parameters, i°f w -0.104I and P^i^f) = 
-0.028FiJ . 

Another important figure of merit is the efficiency, de- 
fined as r\ = — P ou t /-Pin- The plots f?(i ac ) are shown in 
Fig. |U[e). Obviously, the efficiency has a maximum just 
at the beginning of the rectification window and falls with 
increasing |id c |- To derive the approximate behavior of 
r?(J ac ) at the beginning of rectification window analyti- 
cally, we Taylor-expand P ou t and P in near i ac = 7™™(J dc ) 
(|5J). As a result, we get 



T](hc,I 



1 



^dc 



dcJ 



Idc 



6 (Idc " Ic+) 2 



(/ac-i^ (19) 



The first term represents the exact expression for the 
efficiency at the left edge of the rectification window for 
given Idc- To find the maximum efficiency, we vary idc 



from down to IJ^gg 



(I c+ + J c _)/2. At J dc 



the efficiency approaches it ultimate maximum value 

Ic- + Ic+ 



rmax 
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Ic- - I 



c+ 



(20) 



although the rectification window A/ rcct vanishes, see 
Fig. Hljb)-(e). For our values of I c - and I c + (see the 
caption of Fig. HJ) r] max = 0.6. 



V. NON-ADIABATIC DRIVE 

We use a microwave generator with an emitting an- 
tenna close to the ALJJ to drive the ratchet at frequen- 
cies as high as v > 1 GHz. With this geometry, the bias 
leads of the ALJJ act as a pickup antenna. Note that 
the microwaves are also picked up by the injector elec- 
trodes but due to the geometry of the setup it causes only 
a negligible "flashing" part in our ratchet behavior fla ]. 
We average the dc voltage V over 2000 data points at a 
sampling rate of I 00 kHz (~ 10 s periods) and measure it 
vs. the applied power P cx i ac of the generator. 

Rather than growing smoothly as in Fig.^a), the rec- 
tified voltage V is now quantized as 



V, 



sgn (J in ji)n$ ^ 



(21) 



(Shapiro-like steps) |Tl|. Each step corresponds to an in- 
teger number n of turns of a fluxon around the ALJJ per 
one period of the ac drive. The pre-factor — sgn(ii n ji) 
is chosen so that n > corresponds to the motion of a 
fluxon in the "easy" direction of the potential, while n < 
corresponds to the difficult direction (current /voltage 
reversal), cf, Fig. [5J If the time r required for one rev- 
olution of a fluxon around the ALJJ (in the best case 
t~o = Co/ L) becomes comparable with the period of the 
driving force, i.e., if the voltage V = <t>o^ is approaching 
V\, the fluxon has time for only one revolution, so n = I 
and only one step can be observed on the V(P) curves. 
In the non-adiabatic regime, also current reversal can be 
observed, i.e., n = —I resulting in a "negative" voltage 

MM. 



Being in the rectification regime, we apply an addi- 
tional bias current Idc which either tries to stop the 
ratchet operation (i.e., with sign opposite to V) or to 
"help" it (i.e., with the same sign as V). In the non- 
adiabatic regime it makes sense to associate 7 st o P and 
^supp with a particular mode of operation (step n) rather 
than with the ratchet as a whole. Therefore, a par- 
ticular mode of operation is ending when the voltage 
"^Udc ^ 0) ^ V(I dc = 0). The value of I dc at which 
this condition is met is called stopping force i s top (P) or 
support force I sup p (P) respectively. 

Fig. [5] shows a typical rectification curve measured us- 
ing the sample C3 at idc = 0. The voltage step corre- 
sponds to n = — 1. In addition, we show the i s top(P) 
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FIG. 5: Rectification curve with positive potential sign. Stop- 
ping force istop (P) and support force i sup p(P). Black trian- 
gles: P is swept up; grey triangles: P is swept "down". 

and i SU pp (P) dependences for this area. Fig. [5] shows the 
curves for the same parameter setting, only the potential 
amplitude is being reversed in comparison to Fig. [5] 

The rectification curves V(I ac ) show the expected dis- 
crete values of rectified voltage. An inverted potential 
causes an inverted voltage step. Different sweeping di- 
rections of P result in almost identical curves, which sug- 
gests that we are not observing a simple synchronization 
of the fluxon movement with the drive (for which we can 
lock to the n = ±1 step randomly) but rather a true rec- 
tification with well defined direction independent from 
the history. We check this for every voltage step to en- 
sure that we are observing a real ratchet effect and also 
measure i s top and i SU p P for different sweeping directions 
to avoid regimes of synchronization. 

In both plots, i s top = at the lower edge of the rec- 
tification window — a feature that we observe in almost 
all measurements. From this point, the stopping force 
^stop grows smoothly throughout the rectification win- 
dow. It vanishes when we leave the rectification regime 
of n = — 1 . 
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FIG. 6: Rectification curve with negative potential sign. 
Stopping force J s t op (P) and support force I BU pp{P). Black tri- 
angles correspond to sweep up, grey triangles to sweep down. 



The support force / sup p also looks similar in both plots 
(but with opposite sign). It has a minimum value at 
the upper edge of rectification window. At lower driving 
amplitudes / S upp stays ps const, and drops very fast (or 
even jumps abruptly) down to its minimum at the upper 
edge of the rectification window. 

The maximum value of /stop in all our measurements 
was = 215 /iA = 0.112/ (reached at the edge of 

the rectification window) and = 298 fj,A = 0.155/ 

(both in sample C3, not shown). Note, that |/fup P l > 



|/™op I in all measurements. 

In sample E3, we observed both direct rectification 
(n = +1) and reversal (n = — 1) in one rectification curve 
at a driving frequency v = 16.9 GHz, see Fig. Efa). The 
voltage jumps directly from one step to the other, show- 
ing a small hysteresis (indicated by the arrows) when 
sweeping P back and forth. We measured / s top(P) and 
-^supp(-P) for both voltage steps. Also note a data point at 
\J~P ps 0.0048 v"W. The step with n = -1 is metastable 
in the interval \[P = 0.0054 ... 0.0057 \fW, so one can 
see the system jump on it sometimes. 

Fig. EKb) shows I stop (P) and I supp (P) for n = +1. 
The shape of / s to P (P) is similar to the previous plots. 
In the middle of the rectification window, the curve gets 
more noisy and, at the upper edge, / s top shows a small 
hysteresis for different sweep directions of P. I supp (P) is 
almost constant throughout the major part of the plot 
and jumps close to the upper edge of the rectification 
window. 

Fig. EJc) shows /stop and Pupp for the step n = — 1. 
Pupp(P) has a well-known shape — the curve remains at 
high values within the whole rectification window. The 
/stop (P) curve vanishes at both ends of the rectification 
window. In fact, one can see that it consists of two 
branches that join at a value of driving amplitude, where 
the voltage step n = — 1 is ending - i.e., the left branch 
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FIG. 7: (Color online) Rectification curve for v = 16.9 GHz 
(a). Stopping force / st0 p (P) and support force / SU pp (P) for 
voltage step corresponding to n = +1 (b) and n = — 1 (c). 



of /stop (P) shows rectification with n = — 1, although the 
V (P) curve has switched already to n = +1 dynamics 
in this region. The reason is that we are tracking the 
n = — 1 mode, which is not visible for /^c = 0, but ap- 
pears for /dc 7^ 0. The fact that we can track it shows, 
that the presence of the n = — 1 mode depends on the 
history of the system and the value of /d c in a non-trivial 
way. 

The reason for this might be the measurement itself: 
For the measurement of V (P), the current /d c = const, ps 
(a small current offset is always possible) whilst for the 
measurement of / s top (P) or / sup p (P) the current /d c is 
ramped up for each amplitude of the drive P and then 
set to /dc ~ again for the next value of P. This sub- 
tle difference can influence the fluxon dynamics inside 
the junction and therefore the n = — 1 dynamics can be 
pertained longer than what was seen in the rectification 
curve. In this sense, the measurement of the stopping 
force or support force can be used for further investiga- 
tion of possible modes of operation while measuring the 
loading capability of the ratchet. 

For the non-adiabatic regime, numerical simulations 
show a behavior similar to the experiment (for the quasi- 
static regime, these simulations are very time-consuming 
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and can be avoided by using an analytical approach sim- 
ilar to the one discussed in Sec. HVj) . These simulations 
were performed using an explicit numerical scheme for 
Eq. (JTJ) with damping coefficient a — 0.1 (weakly under- 
damped limit). The numerical technique and simulation 
software are discussed in detail in jUJ . 

In Fig. [5] we show simulated (V) (iac) as well as 
^stop(^ac) and Isupp(^ac) dependences. In our simulation 
of V(iac) dependences shown in Fig. [8j the sweep of Z ac 
was performed from 0.3io to 0.6io and back. The value 
of V was calculated by averaging V only over one period 
of ac drive. Therefore, if we are in the chaotic regime, 
we will observe a finite value of the voltage in our sim- 
ulation, while in the experiment these will be averaged 
to some finite value which is often equal to zero. The 
chaotic voltage distribution in Fig. [5] at Z ac > 0.48 shows 
such a chaotic regime. For this simulation, the (discrete) 
step in i ac was chosen very tiny so that the system still 
has time to come in equilibrium, although one sees the 
transient processes on V(I ac ). A large discrete change of 
the system parameters (e.g. the driving power or I^ c ) 
may give a different result than a small one, i.e. the re- 
sulting value of / s top or i SU pp can be much smaller than 
for smooth parameter changes. This happens, because 
the dynamics is close to chaotic and thus even a small 
kick (step-like change of parameters) may drive the sys- 
tem into a different state. Our simulations use a small 
step size of 10 _5 /o- 

Then, for the simulation of i s top(iac) and i SU pp(i a c) we 
started by recalling the state of the system at i ac /io = 
0.425 (where rectification takes place) and then sweeping 
J ac into each direction and, for given J ac , looking for J s to P 
or isupp- For the interpretation of the plots, the curves at 
7 ac < 0.425 have to be interpreted looking at the V(I ac ) 
curve swept in negative direction, while i s top(iac) and 
^supp(-^ac) at / ac > 0.425 have to be followed looking at 
the positive sweep of V(i ac ). 

The V r (/ a c)-curvc reproduces the n = — 1 step sim- 
ilar to the experiment, cf. Fig. [5] and Fig. |6l At 
0.34 < Z ac < 0.41 one can see period two dynamics in the 
V(I&c) curve. Again, this can only be observed in simu- 
lations as these values are averaged to their mean in our 
experiment. One can see that the i s t p(i a c)-curve is sim- 
ilarly shaped as in our measurement (see Fig. [5]) , showing 
an increasing value to larger driving power. At high driv- 
ing amplitudes, one can see a jump in the i s top (iac) curve 
to lower values of i s top which wc did not see in our partic- 
ular measurement. This jump is related to the presence 
of an island in the id c -iac plane (marked by the dashed 
line) where the dynamics is not n = — 1 any more. The 

isupp 

(i ac )-curve shows large values over the whole recti- 
fication window. Wc believe that in our experiment, we 
also meet the island where the solution for n = —1 is not 
stable. 

We have also been able to reproduce a situation in 
which both n = +1 and n = —1 dynamics arc present on 
one rectification curve like in our measurements on sam- 
ple E3. The corresponding numerical simulation results 
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FIG. 8: (a) Numerically simulated for ALJJ of I = 8.0 

and driving frequency ui/uj p \ — 0.25. Black data points cor- 
respond to sweeping / ac upwards, while the gray points to 
a downward sweep. The value of J a c where the simulation 
is restarted (see text) is indicated by the dotted line, (b) 
/stop (Jac) and / sup p(/ac) for the voltage step corresponding to 
n = — 1. Black (gray) symbols correspond to positive (nega- 
tive) sweeping direction of I ac . I s to P is reduced by a region of 
instability (marked by the dashed line). 



for a LJJ of the normalized length I = 10.0 are presented 
in Fig. UK). 

In this simulation, the sweep of I ac was performed also 
from 0.3io to 0.6io and back with the same step in Z ac 
as in Fig. |5Ja). For the interpretation of the plots, for 
Iac < 0.50/rj one has to refer to the V(I ac ) curve swept 
in negative direction, for I ac > 0.50 the positive sweep of 
V(iac) has to be followed. 

In Fig. [9ja) a transition between the voltage steps n = 
+1 and n = —1 show a small hysteresis. Moreover, for 
the negative sweeping direction of I aci a bifurcation to 
"period two" dynamics is visible. In our experiment with 
much longer integration time, this is averaged to the same 
step voltage and is not distinguishable from "period one" 
dynamics. The same applies to the region (outside the 
n = ±1 voltage steps) where the system shows chaotic 
dynamics — due to the averaging in our experiments, 
one observes just an average voltage of the chaotic fluxon 
motion. 

Fig. Ob) shows the I ato p(y/P) and I supP (VP) depen- 
dences for ?i = +1 step. They confirm the behavior ob- 
served in our measurements. Namely, the stopping force 
7 stop (/ac) grows almost linearly from the lower edge of 
the rectification window to larger values and then jumps 
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FIG. 9: (Color online) (a) Numerically simulated V(I ac ) for 
ALJJ of I = 10.0 and driving frequency uj/uj v \ = 0.23. Black 
data points correspond to sweeping / ac upwards, while the 
gray points to a downward sweep, (b) / s to P (/ac) and / S up P (4c) 
for the voltage step corresponding to n — +1. Black (gray) 
symbols correspond to positive (negative) sweeping direction 

Of /ac 



down to zero. The support force / SU pp(^ac) also shows the 
well-known behavior and drops smoothly to smaller val- 
ues while approaching the upper edge of the rectification 
window. 

Fig. Efc) shows the I s to P (-P) and / S u PP (P) dependences 
for n = — 1. They look very similar to those measured 
using the sample E3, c.f. Fig. [TJc). The stopping force 
-^stop (-fac) also grows from the lower edge of the recti- 
fication window and at some point jumps to a branch 
growing from the upper edge - in our simulation, the 
curves also show a jump. This jump is located at a driv- 
ing amplitude / ac where the rectification curve V (/ ac ) 
has a bifurcation and is right inside the hysteretic area. 
The support force / S upp (4c) again drops to lower val- 
ues with increasing driving amplitude I ac and goes down 
to zero at the upper edge of the rectification window. 
One can see that the curve consists of several pieces that 
may be related to the switching onto different dynamic 
regimes. This is in contrast to the experimental results 
where / supp (/ ac ) is falling abruptly down to zero. 

The maximum values I™p = 0.049/n and J™^£ = 
0.167/n. in Fig.^b) and (c) are similar to those obtained 
in our measurements, other simulations also confirm the 



scale of these values. 

Both in simulation and experiment, the maximum 
power delivered to the load in the non-adiabatic regime 
can be obtained by using large driving powers P. As 
shown in Fig. Etc), sometimes the maximum value is 
reached around the center of the rectification window. 
Combining these findings, a good value to start is a large 
value of P on the rectification step with some variations 
for optimization in the actual experiment. Estimated 
maximum P out = / s top^ ~ 50 /j,A x 34 /iV = 1.7 nW for 
n = +1 and P ou t ~ 3.3 nW for n = —1 in Fig. [7] and 
-Pout ~ 5.7 nW for n = —1 in Fig. [6l Note the difference 
from the adiabatic case where V tends to zero smoothly 
as a function of / ac , so that P out — > too. In the non- 
adiabatic case the maximum P out may be reached at the 
edge of the step as in Fig. [51 



VI. CONCLUSION 

We implemented the Josephson vortex ratchet using 
an annular long Josephson junction equipped by current 
injectors used to create the ratchet potential and inject 
a fluxon into the junction. 

In the quasi-static regime we derived the dependence 
of the stopping force 7 s top (lac) analytically using a rea- 
sonable model for the IVC. This model predicts that in 
order to increase the value of the stopping force (and 
other figures of merit), the amplitude of the ratchet po- 
tential should be rather large. The stopping force is di- 
rectly related to the loading capability of a ratchet and 
to the maximum rectified power that it can deliver to the 
load. The maximum output power P out is also reached 
at the upper end of the rectification window for inter- 
mediate values of /<jc- However, the ratchet works with 
maximum efficiency at the lower end of the rectification 
window. The ultimate efficiency that can be reached in 
principle is defined only by the asymmetry of the poten- 
tial, see Eq. (|20|) . 

In the non-adiabatic regime, we have measured recti- 
fication curves and sometimes observed a transition be- 
tween n = +1 and n = — 1 dynamics with a small hystere- 
sis. By applying an additional dc bias current /dc in the 
rectification regime we have measured the stopping force. 
Both / s top and P ut have maximum near the upper edge 
of the rectification window. The support force, i.e., the 
current with the same sign compared to the rectification 
voltage, where the ratchet stops working, is not changing 
very drastically throughout the rectification window. 

Therefore, if one wants to have a robust ratchet op- 
eration tolerant to additional bias currents, the ratchet 
should be operated at the upper end of the rectification 
window and the dc bias current (flowing to the load) 
should not exceed / s top (-^ac )■ 
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